arXiv:1503.04786v5 [math.CA] 1 Jun 2015 


DARBOUX TRANSFORMATIONS FOR 
MULTIVARIATE ORTHOGONAL POLYNOMIALS 

GERARDO ARIZNABARRETA AND MANUEL MANAS 


Abstract. Darboux transformations for polynomial perturbations of a real multivariate measure are 
found. The ID Christoffel formula is extended to the multidimensional realm: multivariate orthogonal 
polynomials are expressed in terms of last quasi-determinants and sample matrices. The coefficients of 
these matrices are the original orthogonal polynomials evaluated at a set of nodes, which is supposed 
to be poised. A discussion for the existence of poised sets is given in terms of algebraic hypersufaces 
in the complex affine space. 


Contents 


1. Introduction 0] 

1.1. Historical background and context [2] 

1.2. Preliminary material 0] 

2. Extending the Christoffel formula to the multivariate realm 0] 

3. Poised sets |9] 

4. Darboux transformations for a general perturbation 06] 

4.1. Discussion for the arbitrary power of a prime polynomial 06] 

4.2. The general case 09] 

References 02] 


1. Introduction 

In a recent paper [H] we studied how the Gauss-Borel or LU factorization of a moment matrix 
allows for a better understanding of the links between multivariate orthogonal polynomials (MVOPR) 
on a multidimensional real space R D , D > 1, and integrable systems of Toda and KP type. In 
particular, it was shown how the LU decomposition allows for a simple construction of the three 
term relation or the Christoffel-Darboux formula. Remarkably, it is also useful for the construction 
of Miwa type expressions in terms of quasi-tau matrices of the MVOPR or the finding of the Darboux 
transformation. Indeed, we presented for the first time Darboux transformations for orthogonal 
polynomials in several variables, that we called elementary, and its iteration, resulting in a Christoffel 
quasi-determinantal type formula. These Darboux transformations allow for the construction of new 
MVOPR, associated with a perturbed measure, from the MVOPR of a given non perturbed measure. 
Observe that they also provide a direct method to construct new solutions of the associated Toda-KP 
type integrable systems. 
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What we called elementary Darboux transformations in pHjj where given as the multiplication of 
the non perturbed measure by a degree one multivariate polynomial. The m-th iteration of these 
so called elementary Darboux transformations leads, therefore, to a perturbation by a multivariate 
polynomial of degree m. This way of proceeding was motivated by the one dimensional situation, 
in that case happens that the irreducible polynomials have degree one (the fundamental theorem of 
algebra). But, in higher dimension the situation is much richer and we do have many irreducible 
polynomials of higher degree. Therefore, the territory explored for the Darboux transformations in 
[H| was only a part, significant but incomplete, of a further more vast territory. In this paper, see 
Theorem 2.1, we give a generalization of the Darboux transformations found in [T4] that holds for a 


perturbation by a polynomial of any degree. This provides us with an elegant quasi-determinantal 
expression for the new MVOPR which is a broad extension of the ID determinantal Christoffel 
formula. 

For the construction of the mentioned general Darboux transformation we use multivariate interpo¬ 
lation theory, see m- Therefore, we need of poised sets for which the sample matrix is not singular. 
In this paper we initiate the study of poised sets for general Darboux transformations. We find that 
the analysis can be splitted into two parts, one measure-independent part depending exclusively on 
the relative positions of nodes in the algebraic hypersurface of the generating polynomial, and an¬ 
other related to the non perturbed measure and the corresponding Jacobi matrices. The geometrical 
part, as usual in interpolation theory, requires of the concourse of Vandermonde matrices. In fact, 
of multivariate Vandermonde matrices, see [42j, or multivariate confluent Vandermonde matrices. 

With the aid of basic facts in algebraic geometry, see for example [3lj or [l6j we are able to show, 
for generating polynomials that can be expressed as the product Q — Q\ ■ ■ ■ Qn of N prime factors, 
-see Theorem 3.3- that there exists, in the complex domain, poised sets of nodes by forbidding its 


belonging to any further algebraic hypersurface, different from the algebraic hypersurface of Q, of 
certain degrees. Moreover, we see that for a perturbation of the measure by a polynomial of the form 



in this general scenario is given in Theorem |4.4[ where again the set of nodes when poised can not 
belong to any further algebraic hypersurfaces of certain type. 

The layout of the paper is as follows. Within this introduction we further perform a number 
observations regarding the historical background and context of the different mathematical issues 
discussed in this paper. Then, we reproduce, for the reader commodity, some necessary material 
from In [|2] we give the Darboux transformation generated by a multivariate polynomial, and 
in ^3] we discuss poised sets, giving several conditions for the nodes in order to constitute a poised 
set. Finally, in [Q we see that the previous construction fails in some cases, and them we present an 
extension of the Darboux transformations which overcomes this problem. 


1.1. Historical background and context. 

1.1.1. Darboux transformations. These transformations were introduced in ([jj5] in the context of the 
Sturm-Liouville theory and since them have been applied in several problems. It was in |4D], a 
paper where explicit solutions of the Toda lattice where found, where this covariant transformation 
was given the name of Darboux. It has been used in the ID realm of orthogonal polynomials quite 
successfully, see for example [I3II91CE1 EM- In Geometry, the theory of transformations of surfaces 
preserving some given properties conforms a classical subject, in the list of such transformations given 
in the classical treatise by Einsehart [22] we find the Levy (Lucien) transformation, which later on 
was named as elementary Darboux transformation and known in the orthogonal polynomials context 
as Christoffel transformation [SHISH]; in this paper we have denoted it by T. The adjoint elementary 
Darboux or adjoint Levy transformation T^ 1 is also relevant [4UJ s20J and is referred some times as a 
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Geronimus transformation [49]. For further information see j45[ 1.30] . For the iteration of elementary 
Darboux transformations let us mention that Szego [4Sj points out that for dp = dx the iteration 
formula is due to Christoffel [TS|. This fact was rediscovered much latter in the Toda context, see 
for example the formula (5.1.11) in [40J for W+{N). 

1.1.2. Multivariate orthogonal polynomials. We refer the reader to monographs |21j and [55]. The 
recurrence relation for orthogonal polynomials in several variables was studied by Xu in [50], while 
in [51] he linked multivariate orthogonal polynomials with a commutative family of self-adjoint 
operators and the spectral theorem was used to show the existence of a three term relation for the 
orthogonal polynomials, ffe discusses in [52j how the three term relation leads to the construction 
of multivariate orthogonal polynomials and cubature formulae. Xu considers in [5_6j polynomial 
subspaces that contain discrete multivariate orthogonal polynomials with respect to the bilinear 
form and shows that the discrete orthogonal polynomials still satisfy a three-term relation and 
that Favard’s theorem holds. The analysis of orthogonal polynomials and cubature formulae on the 
unit ball, the standard simplex, and the unit sphere [53J lead to conclude the strong connection 
of orthogonal structures and cubature formulae for these three regions. The paper pxj] presents a 
systematic study of the common zeros of polynomials in several variables which are related to higher 
dimensional quadrature. Karlin and McGregor [35j and Milch [4TJ discussed interesting examples 
of multivariate Hahn and Krawtchouk polynomials related to growth birth and death processes. 
There have been substantial developments since 1975, for instance, the spectral properties of these 
multivariate Hahn and Krawtchouk polynomials have been studied in [28]. A study of two-variable 
orthogonal polynomials associated with a moment functional satisfying the two-variable analogue 
of the Pearson differential equation and an extension of some of the usual characterizations of the 
classical orthogonal polynomials in one variable was found [ 23] . 

1.1.3. Quasi-determinants. For its construction we may use Schur complements. Besides its name 
observe that the Schur complement was not introduced by Issai Schur but by Emilie Haynsworth in 
1968 in [321133]. In fact, Haynsworth coined that name because the Schur determinant formula given 
in what today is known as Schur lemma in m■ In the book one can find an ample overview on 
the Schur complement and many of its applications. The easiest examples of quasi-determinants are 
Schur complements. In the late 1920 Archibald Richardson |03, B], one of the two responsible of 
Littlewood-Richardson rule, and the famous logician Arend Heyting |34| . founder of intuitionist logic, 
studied possible extensions of the determinant notion to division rings. Heyting defined the designant 
of a matrix with noncommutative entries, which for 2 x 2 matrices was the Schur complement, and 
generalized to larger dimensions by induction. Let us stress that both Richardson’s and Heyting’s 
quasi-determinants were generically rational functions of the matrix coefficients. A definitive impulse 
to the modern theory was given by the Gel’fand’s school |25J [261 EZl EU- Quasi-determinants where 
defined over free division rings and it was early noticed that it was not an analog of the commutative 
determinant but rather of a ratio determinants. A cornerstone for quasi-determinants is the heredity 
principle , quasi-determinants of quasi-determinants are quasi-determinants; there is no analog of 
such a principle for determinants. However, many of the properties of determinants extend to this 
case, see the cited papers. Let us mention that in the early 1990 the Gelf’and school [2E] already 
noticed the role quasi-determinants had for some integrable systems. All this paved the route, using 
the connection with orthogonal polynomials a la Cholesky , to the appearance of quasi-determinants 
in the multivariate orthogonality context. Later, in 2006 Peter Olver applied quasi-determinants to 
multivariate interpolation m, now the blocks have different sizes, and so multiplication of blocks is 
only allowed if they are compatible. In general, the (non-commutative) multiplication makes sense 
if the number of columns and rows of the blocks involved fit well. Moreover, we are only permitted 
to invert diagonal entries that in general makes the minors expansions by columns or rows not 
applicable but allows for other result, like the Sylvester’s theorem, to hold in this wider scenario. 
The last quasi-determinant used in this paper is the one described in [02], see also [ IT] , 
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1.1.4. LU factorization. This technique was the corner stone for Mark Adler and Pierre van Moerbeke 
when in a series of papers where the theory of the 2D Toda hierarchy and what they called the discrete 
KP hierarchy was analyzed CO-0- These papers clearly established -from a group-theoretical setup- 
why standard orthogonality of polynomials and integrability of nonlinear equations of Toda type 
where so close. In fact, the LU factorization of the moment matrix may be understood as the 
Gauss-Borel factorization of the initial condition for the integrable hierarchy. In the Madrid group, 
based on the Gauss-Borel factorization, we have been searching further the deep links between the 
Theory of Orthogonal Polynomials and the Theory of Integrable Systems. In |8j we studied the 
generalized orthogonal polynomials [lj and its matrix extensions from the Gauss-Borel view point. 
In pi] we gave a complete study in terms of factorization for multiple orthogonal polynomials of 
mixed type and characterized the integrable systems associated to them. Then, we studied Laurent 
orthogonal polynomials in the unit circle trough the CMV approach in |TUJ and find in m the 
Christoffel-Darboux formula for generalized orthogonal matrix polynomials. These methods where 
further extended, for example we gave an alternative Christoffel-Darboux formula for mixed multiple 
orthogonal polynomials H2 or developed the corresponding theory of matrix Laurent orthogonal 
polynomials in the unit circle and its associated Toda type hierarchy [13]. 

1.2. Preliminary material. Here we remind the reader some necessary content extracted from 
m Our method to construct Darboux transformations of multivariate orthogonal polynomials in 
a D-dimensional real space (MVOPR) is formulated in terms of a Cholesky factorization of a semi¬ 
infinite moment matrix. We consider D independent real variables x = (xi,x 2 , . ..,Xd) E 0 C M D 
varying in the domain together with a Borel measure d fi(x) E The inner product of two 

real valued functions f(x) and g(x) is defined by 

if,g) '■= [ f( x )g( x )dfi{x). 

Jo. 

Given a multi-index at = (aq,..., ctn) T E HJff of non-negative integers we write x a = x " 1 • • • xff } ; 
the length of a is |a| := a a- This length induces the total ordering of monomials, x a < x a ' 

| o: | < | a'|, that we will use to arrange the monomials. For each non-negative integer k E Z + we 
introduce the set 


[k] := {a E Zif : |a| = k}, 

built up with those vectors in the lattice Z+ with a given length k. We will use the graded lexico¬ 
graphic order; i.e., for aq,CK2 E [k] 


aq > aq 4=> 3p E Z + with p < D such that a u = a 2 ,i, ■ ■ ■, aq, p = a 2jP and aq jP+ i < a 2 , P + 1 , 

and if at.^ E [k] and E [^], with k < i then < cPA Given the set of integer vectors of 

length k we use the lexicographic order and write 

[k\ = {af } , a4 fc) , • • •, aj^} with > a%. 

Here | [k] \ is the cardinality of the set [k], i.e., the number of elements in the set. This is the dimension 
of the linear space of homogenous multivariate polynomials of total degree k. Either counting weak 
compositions or multisets one obtains \{k}\ = ((f)) = ( D+ ^~ 1 )- The dimension of the linear space 
Mfc[xi,..., xjj] of multivariate polynomials of degree less or equal to k is 

N k — 1 + I [2] | "I-f | [k] | = ^ D ^ • 
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We introduce the vector of monomials 

/X[o]\ 

X[i] 


X ■ = 


X[k] 

\ : / 


/ x ai \ 


where 


X[k] 


x 


«2 


\a; a l[ fe ]l ) 


Observe that for k — 1 we have that the vectors a A = e a for a G {1,..., D} forms the canonical 
basis of M d , and for any cxj G [k\ we have olj = ^ 0=1 a“e a . For the sake of simplicity unless needed 

we will drop off the super-index and write ctj instead of ot^ k \ as it is understood that |a,-| = k. 

The dual space of the symmetric tensor powers is isomorphic to the set of symmetric multilinear 

functionals on M D , (Sym fc (M D )) = S'((M D ) fc ,M). Hence, homogeneous polynomials of a given total 

degree can be identified with symmetric tensor powers. Each multi-index a G [k] can be thought as a 

weak ^-composition of k (or weak composition in D parts), k = a± + - ■ - + od. Notice that these weak 

compositions may be considered as multisets and that, given a linear basis {e a }^ =1 of we have the 

linear basis {e ai 0- • •©e afc }i<ai<---<a t .<T> for the symmetric power S fc (M D ), where we are using multisets 

fcez + 

1 < op < • • • < a*, < D. In particular the vectors of this basis e ® M (“d q ... © 6 © m (°p), or better its 
duals 0 • • • 0 (e*J GM( “ p ) are in bijection with monomials of the form • • ■ x^} ap \ 


The lexicographic order can be applied to (M d ) G/c = we then take a linear basis of S fc (M D ) as 


the ordered set B c = {e ai ,..., e“iwi} with := © • • • 0 e^ 3 so that X[k\(x) = X^=i x aj e aj . 

We consider semi-infinite matrices A with a block or partitioned structure induced by the graded 
reversed lexicographic order 


OcA 



( A io\m 

Ao],[i] • • A 



• * A oc ^ ^ 

A= j 

A mo] 

K ■ 

di[l],[l] • • • 

: ) 

| » Akim = 

\ a |[fc]|’ a l 

. . A (fc) (£) , 
a m\> a \[i\\/ 


We use the notation o [fc]lM e mm\m for the rectangular zero matrix, Orw G for the zero vector, 
and I [fc] G for the identity matrix. For the sake of simplicity we normally just write 0 or I 

for the zero or identity matrices, and we implicitly assume that the sizes of these matrices are the 
ones indicated by its position in the partitioned matrix. 

Definition 1.1. Associated with the measure d n we have the following moment matrix 


G: = X(*)d h(x) X (x) 


T 


We write the moment matrix in block form 

' G [ ;o],[< 

G = I G [i],[i 


G [o],[i] 

G [i],[i] 


Truncated moment matrices are given by 


G m : = 


G u 


\G\f- 


[<-!],[' 
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Notice that from the above definition we know that the moment matrix is a symmetric matrix, 
G = G T , which implies that a Gauss-Borel factorization of it, in terms of lower unitriangular [] and 
upper triangular matrices, is a Cholesky factorization. 


Proposition 1.1. If the last quasi-determinants 0*(G[ fc+1 l), k e {0,1,. 
matrices are invertible the Cholesky factorization 

(1.1) G = S~ 1 H(S~ 1 ) T , 

with 


/I 0 0 ■ ■ A 

OS" 1 ) HUP] I 0 ••• 

( S 1 )[2],[0] (S 1 )[2],[1] I 

v •••/ 


/ 77 101 
0 
0 

V i 


} , of the truncated moment 


0 0 \ 
H m 0 

0 H[ 2 ] 

•7 


can be performed. Moreover, the rectangular blocks can be expressed in terms of last quasi-determinants 
of truncations of the moment matrix 


H [k] = e*(G[ fc+1 i), (s-'hm = e*(G ! if +1] )e*(G [m] )- 1 . 


We are ready to introduce the MVOPR 


Definition 1.2. The MVOPR associated to the measure dp, are 
fRr.A , ( P JH>\ 


(1.2) p = Sx= P hl , P[k]{x) = y^5[fci,[/|X[<|(g) 


£=0 


\ a l[fe]l/ 


• trES s , 


1=0 j=1 


(£) 

(fc) {t)X a i . 

ot ; Got) 

i » j 


Observe that P [ k ] = X[k]( x ) + P[k]X[k-i]( x ) + ■■■ is a vector constructed with the polynomials 
P a .{x) of degree k , each of which has only one monomial of degree k ; i. e., we can write P az (x) = 
x a% + Q ai ( x ), with deg Q ai < k. 


Proposition 1.2. The MVOPR satisfy 


(1.3) 

(1.4) 


P[k]( x )dp(x)(P [£] (x)) T = I P [k] (x)dp{x)( X [£](x))' =0, 


= 0,1, • ■ ■, k - 1, 


in 


P [k] ( x )dp(x)(P [k] (x)) T = I P [k] (x)dp(x)(x[k]( x ))' = H [k] . 


in 


Therefore, we have the following orthogonality conditions 




/ P a w( x )P a w( x ) dp(x)= / P <k) (x)x a i dp(x) = 0, 

Jn i j Jn i 

for t = 0,1,..., k — 1, i — 1,..., | [k] and j = 1,..., | [I] |, with the normalization conditions 

[ P ai ( x )P aj ( x ) dp(x)= [ P a .(x)x aj dp(x) = H a . !a ., i, j = 1,..., \[k]\. 

Jn Jn 

Definition 1.3. The shift matrices are given by 

/0 (A a ) [0 ],[i] 0 0 ■■■\ 

0 0 (A a )[i],[ 2 ] 0 

00 0 (A a ) [ 2 ],[3] 

0 0 0 0 


A a — 


\ = 


7 


^Lower triangular with the block diagonal populated by identity matrices. 
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where the entries in the non zero blocks are given by 

(-^a) ^a.^ +e a ,a.( k+1 ' 1 ’ ® !)•••) ^i ^ 1) ■ ■ ■ i I [^] 1 i 3 1, . . . , | [fc “I - l] |, 

and the associated vector 

A := (Ai,..., Ad) t . 

Finally, we introduce the Jacobi matrices 

(1.5) J a :=SA a S~\ ae{l,...,D}, 

and the vector 


J = (Ji,..., J d ) t . 

Proposition 1.3. (1) The shift matrices commute among them 

A a A b = A b A a . 

(2) We also have the spectral properties 

(1.6) A a x(x) = x a x(x). 

(3) The moment matrix G satisfies 

(1.7) A a G = G(A a ) T . 

(4) The Jacobi matrices J a are block tridiagonal and satisfy 

J a H =HjJ , ae{l,...,D}. 

Using these properties one derives three term relations or Christoffel-Darboux formulae, but as 
this is not the subject of this paper we refer the interested reader to our paper [14]. 

2. Extending the Christoffel formula to the multivariate realm 

In this section a Darboux transformation for MVOPR is found. Here we use polynomial perturba¬ 
tion of the measure, but to ensure that the procedure works we need perturbations that factor out 
as N different prime polynomials. Latter we will discuss how we can modify this to include the most 
general polynomial perturbation. 

Definition 2.1. Given a degree m polynomial Q e Mfcc], deg Q = m, the corresponding Darboux 
transformation of the measure is the following perturbed measure 

d n(x) Td n(x) := Q(x) d p,(x). 

Observe that, if we want a positive definite perturbed measure, we must request to Q to be positive 
definite in the support of the original measure. 

From hereon we assume that both measures d /i and Q d /i give rise to well defined families of 
MVOPR (equivalently that all their moment matrix block minors are nonzero). 

Proposition 2.1. IfTG is the moment matrix of Tdp we have 

TG = Q(A)G = G{Q{ A)) T 

Proof. It is a direct consequence of the spectral property Q(A)x{x) = Q(x)x(x). □ 
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Definition 2.2. The resolvent matrix is 

to :=(TS)Q(A)S~ 1 

given in terms of the lower unitriangular matrices S and TS of the Cholesky factorizations of the 
moment matrices G = S' -1 i/(S' -1 ) T and TG = (TS)~ 1 (TH)(TS~ 1 Y. The adjoint resolvent is 
defined by 

M := S(TS)-\ 


Proposition 2.2. We have that 

( 2 . 1 ) 


cTH)M t H = u 


Proof. It follows from the Cholesky factorization of G and TG and from (2.1). 

Proposition 2.3. In terms of block superdiagonals the resolvent co can be expressed as follows 

co= Q (m \ A) 

m-th superdiagonal 

+ (T/3)Q (m-1) (A) - Q (m - 1} (A)/3 

(m — 1 )-th superdiagonal 

+ ( TH)H ~ 1 


□ 


diagonal 

Proof. The adjoint resolvent is a block lower unitriangular and the resolvent lo has all its superdiag¬ 


onals but for the first m equal to zero. The result follows from (2.1). 
Proposition 2.4. The following LU and UL factorizations 

Q{J) =Mco, Q{TJ ) =uM, 

hold. 


□ 


Proof. Both follow from Proposition |2.1| and the Cholesky factorization which imply 

(TS)~ 1 H(TS- 1 ) t = Q(A)S~ 1 H(S~ 1 ) T , 

and a proper cleaning do the job. □ 

From the first equation in the previous Proposition we get 
Proposition 2.5. The block truncations (Q(J))^ admit a LU factorization 

(Q(J)) [k] = M [k] cu [k] 

in terms of the corresponding truncations of the adjoint resolvent and resolvent uf k ^. 


Proposition 2.6. We have 


det(Q(J)) [fc] 


and therefore (<2(J))^ is a regular matrix. 


^ det TH\i] 

In det ^l 


Proof. To prove this result just use Propositions |2.4| and |2.3| and the assumption that the minors of 
the moment matrix and the perturbed moment matrix are not zero. □ 


Proposition 2.7. The MVOPR satisfy Q{x)TP{x ) = uiP(x). Consequently, for any element p in 
the algebraic hypersurface Z(Q ) := {x 6 : Q(x) = 0} we have the important relation 

(2.2) ^[fc] \k+m}P[k+m\ i,P) T ^[k],[k-\-m— l]P[k-\-m— 1] (p) d - ’ ’ ’ T ^[k],[k]-P[k] ip) 0- 
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Proof. We have 

uP(x) = (: TS)Q(A)S~ 1 S X (x ) 

= {TS)Q( A) X (x) 

= Q(x)(TS) X (x) 

= Q(x)(TP)(x). 

Finally, when this formula is evaluated at a point in the algebraic hypersurface of Q we obtain 
that the MVOPR at such points are vectors in the kernel of the resolvent. □ 

To deal with this equation we consider 

Definition 2.3. A set of nodes 

M k ,m ■= {VjY-=T C 

is a set with rk, m = Nk+m-i ~ Nk-i = | [&] | 4— • + | [k + m — 1] | vectors in R D . Given these nodes we 
consider the corresponding sample matrices 

( P [k](Pl) ••• P [k](Pr k ,J \ 




e M rfe ’ 


m ^ Pk.m 


yPlfc+TO—1] (Pi) ... P[k+m— 1] (Pr|., m )/ 

S[»,H : =UWPi)> • • •. P[*+™)(Pn.J) e RlI ‘’ +m " 


xr fe . 


Lemma 2.1. When the set of nodes N k ,m C Z(Q) belongs to the algebraic hypersurface of the 
polynomial Q the resolvent coefficients satisfy 


k],[k-\-m] Pj [k,m] T (^[fc],[fc]; ■ ■ ■ j ^[k],[k+m— O’ 


□ 


Proof. Is a direct consequence of (2.2). 

Definition 2.4. We say that Mk, m is a poised set if the sample matrix is non singular 

det E™ ± 0. 

Theorem 2.1. For a poised set of nodes Mk, m C Z(Q) in the algebraic hypersurface of the generating 
polynomial Q the Darboux transformation of the orthogonal polynomials can be expressed in terms of 
the original ones as the following last quasi-determinantal expression 


TP ' k](x) = - 0 ( 5 )- 




0 * 




V S 


[k,m\ 


P[k](x) \ 

P[k+m— 1] (*e) 
P[k+m] (®) ) 


Proof. Observe that Lemma 2.1 together with uJ[k],[k+m} = (Q(A)) 

{U[k],[k], ■ ■ ■ ,U>[ k \ t [k+m- 1]) = — (2(-^-))[fc] j [fc +m 

and from Q(x)TP(x ) = ooP(x) the result follows. 


[&],[&+ 


implies 

-i 


□ 


3. Poised sets 

To construct Darboux transformations in the multivariate setting we need of poised sets in order 
to find invertible sample matrices with the original polynomials as interpolating functions. When is 
this possible? Let us start a discussion on this question. First, we introduce two important matrices 
in the study of poised sets 
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Definition 3.1. We consider the Vandermonde type matrix 

V? ■= (x [i+m HpX .... x |t+m, (p,,,J) £ H"**”*-' 


made up of truncated of multivariate monomials ~xf k+m \x) evaluated at the nodes. We also consider 
the following truncation S™ G W' k ’ mXNk+m - 1 of the lower unitriangular factor S of the Gauss-Borel 
factorization of the moment matrix 


(3.1) 


nm . _ 

°k — 


/ s mo] 

%+l],[0] 


S[k],[i] ■ ■ • %]| 0[fc],[fc+i] 

<%+!],[i] • • • %+i],[fc] %+i]| 


0[fc],[fc+m— 1] \ 
0[fc+l],[fc+m—1] I 


1], [0] $[k-\-m— 1], [1] 


s 


[k+m— 1], [k+m— 2] 


I|[fc+m—1]| 


/ 


They are relevant because 


Lemma 3.1. We have the following factorization 

Y ^ 771 _ Qmym 

From where it immediately follows that 


Proposition 3.1. The following relations between linear subspaces 

Ker Vr C Ker E™, Im E™ C Im SJf = IT, 

hold true. 


The poisedness of Af k , m can be reworded as 

KerEr = {0}, 

or equivantlently 

dim Im E™ = 

Proposition 3.2. For poised set A4, m the multivariate Vandermonde V™ needs to have full column 
rank; i.e., dim Im VJf' = r^ m . 

Proof. For a set to be poised we need that KerS™ = {0}, but KerV™ C KerE™ and consequently, 

dim Ker V™ = 0, 

and, as dim Ker VJf + dim Im VJf' = rk, m , full column rank of the Vandermonde matrix is needed for 
a set to be poised. □ 

The study of the orthogonal complement of the rank; i.e, the linear subspace (Im V™) ± C 
of vectors orthogonal to the image Im V™ where v G (imV™) 1 if v T V™ = 0, gives a better insight 

on the structure of the rank of the Vandermonde matrix. As Im VJf' © (Im = M Ar * : + m ' 1 we have 

the dimensional formula 

dim (Im V™) -1 + dim (Im VJf 1 ) = N k+m _ 1 . 

Proposition 3.3. The Vandermonde matrix Wff has full column rank if and only if 

dim (Im VJf ) ± = N k+m _i - r Km = N k _ i. 

Algebraic geometry will guide us in the search of poised sets in the algebraic hypersurface of 
the m-th degree polynomial Q, Af k , m © Z(Q). We need to abandon the real field M and work in 
its algebraical closure C; i.e., we understand Q as a complex polynomial with real coefficients and 
consider its zero set as an algebraic hypersurface in the D-dimensional complex affine space C D ; we 
also change the notation from x G M. D to z E C D . 
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Definition 3.2. For a multivariate polynomial V of total degree deg V < k + m — 1 its the principal 
ideal is (V) := C {z a V(z) : a G Z+} C C[z±, ..., zjf\ and for its intersection with the polynomials of 
degree less or equal than k + m — 1 we employ the notation 

(V^) k+m— 1 = (v)nc k-\-m— 1 [zi, ...,z D } = C{z a V(z)}o<\ a \ <k+m —deg V • 

It happens that the elements in the orthogonal complement of the rank of the Vandermonde matrix 
are polynomials with zeroes at the nodes 

Proposition 3.4. As linear spaces the orthogonal complement of the rank of the Vandermonde matrix 
(Im VJf j and the space of polynomials of degree less than k + m and zeroes at Afk, m are isomorphic. 

Proof. The linear bijection is 

Nk+m—1 1 

v = K)£ 0 +m " 1 “ 1 G (Im V™) ~ L GG V(z) = v * z<Xi 

i=0 

where V(z) does have zeroes at Af™ C Z(V). Now, we observe that a vector v = 1 G 

(Im VJf ) can be identified with the polynomial V ( z) = JV ViZ OLi which cancels, as a consequence 
of n T V™ = 0, at the nodes. □ 


Thus, given this linear isomorphism, for any polynomial V with deg V < k + m with zeroes at A f™ 
we write V G (imVf) 1 . 

Proposition 3.5. Given a polynomial V G (Im V™) X then 

Mto-i c (Imvr)f 

or equivalently 

(V)+-, 2 Imvr- 


Proof. Observe that all the polynomials in {V)k+ m -i have zeroes at the nodes and have degree less 
or equal to k + m — 1; hence, we have for the corresponding vectors {v a } C (Im V™) . □ 


An important result is 

Theorem 3.1. Given a polynomial Q, deg —m, and a set of nodes Mk, m the corresponding Vander¬ 
monde matrix V™ has full column rank if and only if 

(swi = (imtvny. 

Proof. We know that dim(Q)/ c+m _i = Nk-i and as A4, m C Q we know that (Q)^ +m _ 1 5 Irn VJf'. 
Thus, full column rankness of the Vandermonde matrix can be achieved if only if (Q)k+m- 1 = 

(ImW)) 1 . □ 

From the spectral property Q(A)x{z) = Q(z)x(z) we deduce 

Proposition 3.6. The row (Q(A)) q ,, a G Z^ 7 is the graded lexicographic ordering of the entries in 
the corresponding polynomial z a Q(z). 

Thus, in some way Q( A) encodes the same information as the principal ideal of Q does. To make 
this observation formal we first consider the 
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Definition 3.3. The matrix (Q(A))[ fe,m l G C Nk - lXrk ’ m is given by 


(Q(A))[ fc,m l := 


^ (Q( A ))[o],[fc] 

■ (2( A ))[0],H 0[0],[m+l] 

0[0],[fc+m—1] 

(Q( A ))[fc-2],[jfc] 


0[fc—2],[/c+m— ‘ 

\(Q(A))[ fc -i ]l[fc] 


••• (Q( A ))[fe-l],[fc+r 

/ °[0],[fc] 

°[0],[fc+l] 

0[0],[fc+m—1] ^ 

^[k—m—l],[k] 

®[k—m— l],[fc+l] • • • 

0 [k—m— 1], [k-\-m— 1] 

( Q(A) )[fc_m] ,[fc] 

m],[/c+l] • • • 

0 [k—m ], [k+m— 1] 

(Q(A))[fc_ m _|_i] [j.] 

( Q ( A ) ) [fe—m+l], [fc+1] 




0 [k— 2], [k-\-m— 1] 

\ (Q(A))[ A ._i] j [ fc ] 

(Q( A ))[fc-i],[fc+i] 

(Q( A ))[fc-l],[fc+m-l]/ 


We collect this matrix and the truncation (Q(A))^ in the N^-i x Nk+ m -i matrix 

(G(A))? := ((C(A))I*1, (C(A))I‘'”1). 


k < m, 


k > m. 


Proposition 3.7. We have the following isomorphism 

P= PcxZ a e (Q)k+m-l (-Pao, ■ ■ ■ , = («(), ■ ■ ■ , Ojv # ._ 1 )(Q(A))^ 1 . 

\cx\<k+m 

between the truncated ideal (Q)k+m-i an d the orbit of C Nk -' under the linear morphism (Q(A))™. 
Here we have ordered the multi-indices a in P according to the graded lexicographic order ck 0 < • • • < 

k+rn -l ■ 

Proof. Just recall that P is going to be a linear combination of the polynomials z a Q(z), |a| < k. □ 

Now we show that full column rankness of the Vandermonde matrix is not only a necessary but 
also a sufficient condition 


Theorem 3.2. Let Q G C[z\ be a m-th degree polynomial, deg Q — m. Then, the set of nodes 
Mk,m C Z(Q) is poised if and only if the Vandermonde matrix V™ has full column rank. 

Proof. Let as assume the contrary, then the sample matrix E™ is singular, and we can find a nontrivial 
linear dependence among its rows (E™)j, i G {1,... , r fc m } of the form 

r k,m 

I>(s?), = o, 

i=l 


for some nontrivial scalars {ci,..., c rk m }. But, according to Lemma 3.1 (E™)j = (Sjf)iV™, where 
(S™)i is the i-th row of 5^ and we can write 

r k,m 

E c i(A ra ).W = o, 


2—1 


so that 


rk,-. 


J2^sr),e(imvrr, 


i— 1 


and given the column full rankness of the Vandermonde matrix, see Theorem 3.1, we can write 

(ci,...,c r ,,jsre(imvr)h 
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or, following Proposition |3.7| we get 

(3-2) (a 0 , • • •, a J v fc _ 1 )«2(A))3T = (<*, ■ • •, 

for some non trivial set of c’s. This was the first part of the proof. 
For the second part, we focus on the relation 

SQ( A) = Q(J)S. 

Recalling the S is lower unitriangular by blocks 


S = 


SM I 0 

0 ... 

ST 

0 ■■■ 

* 


and that either Q(A) or Q(J) are block banded matrices with only m block superdiagonals non zero 


2(A) = 


(Q((A))W (Q((A))[ fc ’ m l 0 


Q(J) = 


(2((J)) [fc] (2((J)) [fc ’ m l 0 


* * > 
which can be written as follows 

(3.3) S [fc] (Q(A))™ = (Q(J)) [fe ’ m] 5r+((2(J)) [fc] ^ 1 ,0). 

We now assume that ( |3.2| ) holds and multiply ( |3.3 ) by its left with the nonzero row vector 

O^.KSW )" 1 

to get 

K ■ ■ ■, o Wl _J(e(A))” = (a* ... .o f , 1 _,)(sW)- 1 (a(J)) l ‘ H 5r 

+ (%.% i .,)(s l ‘ , )' 1 ((S(J)) l ‘ l s [ ‘ l .o). 


The expression (3.1) can be written as S™ = (Sj! 1 , S ™) where first block S™ is a r^ m x N k _ x matrix 
and the second block S™ is a lower unitriangular r k ,m x rk, m matrix. Hence, the system (3.2) can be 
splitted as follows 

(«o,... ,o Wl .,)(SW)- I (a(J)) M Sr + ( 00 ,.... 0 W ,. 1 )(5l*l)- 1 (i2(J)) l * l (S w )" 1 = (c„.. .,c ri JST, 

(oo, ■ • • ,o w ,_ 1 )(SW)- I (2(J)) [ ‘'"‘ 1 Sr = (ci, • • • ,c n JS?. 

As S™ is unitriangular it is invertible and from the second equation we get 

(oc-.-.o^.J^Wj-'OtJ)) 1 *'” 11 = (C, .c ri J. 


Thus, (3.2) holds if and only if 


(oo,...,o Jfl _,)(SW)" I (C(J)) [ ‘ 1 SW = o. 


or equivalently if and only if det(Q(J))^ = 0. Thus, recalling Proposition 2.6 (and our inital 
assumption that the two measures d fi(x) and its perturbation Q d /i{x) do have MVOPR) the results 
follows. □ 


Theorem 3.3. Let Q = Qi ■ ■ ■ Qn be the product of N different irreducible polynomials with deg Q a = 

N N 

m a , a G {1,..., N}, and deg Q — m— m a■ Then, the set A4,m C Z(Q) = |J Z(Q a ) is poised if 

a =1 a =1 

the nodes dot not belong to any further complex algebraic hypersurface of degree smaller than k + m 
and different from Z(Q). 



















14 


GERARDO ARIZNABARRETA AND MANUEL MANAS 


Proof. Given a subset Y c C D we define the corresponding ideal I(Y) = {P G C[zi,..., zjf\ : P{z) = 

N N N 

OVz G Y}\ then, /( IJ Y a ) = fl ? (W) and therefore I{Z{Q )) = f) I{Z(Q a )). But, according to 

a= 1 a= 1 a= 1 

the Hilbert’s Nullstcllensatz and the prime character of each factor Q a (every prime ideal is radical) 
we can write 

N 

me))=V®=r\m 

a= 1 

=« 2 ) 

where y/ ( Q ) is the radical of the principal ideal of Q. Thus, we conclude 

(Im V ™)" L D (Q) k+m - 1 

and deduce 

dim (imV™) 1 > N k ~i. 

The equality is achieved whenever we can ensure that there is no further algebraic hypersurface 
of degree less than k + m, different from Z(Q), to which the nodes also belong; i.e., lm(V™) = 

(Qi ■ ■ ■ n 

We now discuss on the distribution of nodes along the different irreducible components of the 
algebraic hypersurface of Q. 

Proposition 3.8. In a poised set, D > 1, the number n a of nodes in the irreducible algebraic 
hypersurface Z{Q a ) fulfill 

k T m a Y n a Y r k+m—m a ,m a i ^1 T ■ ■ ■ T ujv f k,m■ 

Proof. Assume that a number M a smaller than r a k + m a , M a < r a , of nodes lay in the ir¬ 
reducible algebraic hypersurface Z(Q a ), i.e., the number of nodes in its complementary algebraic 
hypersurface Z(Qi ■ ■ ■ Q a _iQ a+1 ■ ■ ■ Q N ) is bigger than r km — r a . Then, the set of nodes belong to 
the algebraic hypersurface -different of Z(Q)~ of degree m — m a + M a < k + m of the polynomial 
Q i ■ • ■ Q a _iQ a+ i • • ■ <2 at7Ti • • • m a , where ttj is a degree one polynomial with a zero at the j'-tli node 
that belongs to Z(Q a ), where we have taken care that Ti\ ■ ■ -7 tm u ^ {Qa), which for D > 1 can be 
always be achieved. Therefore, we need M a > k + m a to avoid this situation and to have a poised 
set. 

The maximum rank of the Vandermonde submatrix built up with the columns corresponding to 
the evaluation of x a t the nodes in Z(Q a ), recalling that dim(Q a )/ c+m _i = Nk+ m -m a , is N k+rn -\ — 

Affc+m— 1— m a fk+m—m a ,m a ■ I—I 

Notice, that we need to put k + mt nodes at each irreducible component Z(Q a ), for a G {1,..., N}, 
hence we impose conditions on Nk + m nodes. But, do we have enough nodes? The positive answer 
for D > 1 can be deduced as follows. 

Proposition 3.9. The bound r k m > Nk + m holds. 

Proof. We have r k , m — \[k + m — 1]| -I-+ |[fc]|, thus a rude lower bound of nodes (for D > 1) is 

r k ,m >m\[k]\ =m(^ + D D _ ^ ^ = m(l + ■ ■ • (l + (1 + k) 

>m(k + 1) 

>Nk + m. 

□ 
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But, what happens with this condition for D — 11 Now, we have r k , m — m nodes, m a = 1 , each 
Z(Qi) is a single point in C and N = m. In this case, the reasoning that lead to the construction of 
the polynomial Qi ■ ■ ■ Q a -iQ a +i ■ ■ ■ Qn^i • • • in the previous proof is not applicable; first M a = 1, 
given that all prime factors are degree one polynomials and second, the polynomial 7Ti must be Q a 
and therefore the product leads to the polynomial Q and no further constraint must be considered. 

The maximum n a is greater than the minimum number of nodes of that type rk+m-m a ,m a > 
m a {k + m — m a + 1) > k + m a . Moreover, the sum of the maximum ranks exceeds the number of 
nodes, and the full column rank condition is reachable: 

N 

Proposition 3.10. We have YL r k+m-m a ,m a > r k,m■ 

i= 1 

Proof. For iV = 2 we need to show that r k+m . 2 , mi + r k+miim2 > r k)Tni+m2 or 


| [k + m 2 + mi - 1] | 4- b | [k + m 2 \ \ + \ [k + m 2 + mi - 1] | 4-h | [k + mi] | 

> | [k + m 2 + mi - 1] | 4-h | [k + m 2 \ \ + \ [k + m 2 - 1]\ -b | [k] 


which is obvious. Then, for N = 3 we need to prove that r k+m2+m3 , mi +r k+rni+rn3 . m2 +r k+mi+m2tm3 > 
h,mi+m 2 +m 3 j but using the already proven N = 2 case we have r k+m2+m3m + r fc+mi+m3 , m2 > 
r k+ m 3 ,m 1 +m 2 and using the N = 2 equation again we do have r k+m3;m , +m2 +r k+m , +m2)m3 > r k)mi+m2+m31 
as desired. An induction procedure gives the result for arbitrary N. □ 


In the next picture we illustrate the case N — 2 of two prime polynomials of degrees rni and m 2 . 
The blue rectangle gives the possible values for the number of nodes ( 711 , 712 ) corresponding n l to the 
prime polynomial Q t according to the bounds 


k ~b mi ^ Tl\ f T' k +m, 2 .m,\ 7 


k "b m 2 A 71 2 A fk+m\,m2- 


The blue diagonals rq + n 2 = K are ordered according (we assume for the degrees that mi < m 2 
and therefore r k+m2jmi < r k+rrl]:m2 ) to the chain of inequalities 


2 k ~b m ^ max(rjfc +m2>mi ~b k T m 2 , iT+mi,m 2 ~b A: ~b mi) A i" k ,m — h+m 2 ,mi T iVc+mi .m 2 • 


Notice that rnax(r fe+m2)mi + k + m 2 ,r fc+mijm2 + k + mi) < r fe>m follows r fcjm = r k+mumj + > 

Tk+mum, + 7Uj|[A:]| > r k+mumj + m^/c + 1) > r k+mum . + k + TTi.j, where (i, j) = (1, 2), (2,1). Therefore, 
the striped triangle is the area where the couples ( 711 , n 2 ) of number of nodes belong. We have drawn 
the passing of the line n\ + n 2 = m trough it, and show the integer couples in that segment, those 
will be the possible distributions of nodes among the zeroes of both prime polynomials. 
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1 %i — k -\~ 7Tli Til 



4. DARBOUX TRANSFORMATIONS FOR A GENERAL PERTURBATION 


We begin with a negative result 

Proposition 4.1. Poised sets do not exist for Q = TZ d , d G {2,3,... } , for any given polynomial 1Z. 

Proof. Now Q = TZ d , d G {2,3,4,...} deg Q = d deg 77, for some polynomial 77. In this case 
Z(Q) = Z(TZ), but dim(7^)fc +m _i = Nk-i+(d-i)degn > -W and consequently the set is not poised. □ 

We now discuss a method to overcome this situation. We will generalize the construction of nodes, 
sample matrices and poised sets. In this manner we are able to give explicit Christoffel type formulae 
for the Darboux transformation of more general generating polynomials. We consider multi-'Wronski 
type matrices and multivariate confluent Vandermonde matrices. 


4.1. Discussion for the arbitrary power of a prime polynomial. Now we take Q = lZ d , 
deg 1Z = n and deg Q = dn , so that Z(Q) = Z{1Z) with 1Z to be a prime polynomial. From 
Proposition 2.7 we know that uP(x) = TZ cl (x)TP(x). To analyze this situation we consider a set of 


linearly independent vectors {nf j P d_ 


C 


D ) &J , Pj < | [j] |; here n [ f ] = (n^) a6 p], and to 


O') _ („ 0) ■ 


each of these vectors we associate the following homeogenous linear differential operator 


d l = V n (i) P_ 
dnP ,V 

* \<X\=] 


From the Leibniz rule we infer 
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Proposition 4.2. For any element p in the algebraic hypersurface Z(TZ) : = {x G R D : TZ(x) = 0} 
we have 


d j P r 


w ^ [K-\-na\ / \ 

^[k\,[k+nd] Tj -j vJP/ ' 1] 

for j G {0,1, ..., d — 1} and i e Pj}. 


(FP- 


[k+nd— 1] 


d j P< 


ip) + ■ ■ ■ + = °> 


dn 


00 


dn 


00 


This suggests to extend the set of nodes and the sample matrices 
Definition 4.1. We consider the splitting into positive integers r^nd — Nk+nd-i — N k _± = \[k] \ + 

d -1 Pj ... 

• • • + | [k + dn — 1] | = h i an d f or eac h j £ {0,1,..., d — 1} we consider the following set of 

j =0 i =1 

distinct nodes 


(j) 

K U) : = {PmKIi c 




where we allow for non empty intersections between these sets of nodes and we denote its union 

d -1 Pj n 

by A 4, U U-A/T. the also need of the above mentioned set of linearly independent vectors 

3=0 *=1 

C = (IRO 0 ) 0 ^ pj < | [j] |, j G {0,..., d — 1}. The partial blocks of the homogeneous 
sample matrices are 


/ ^ (pg ; 




nd\(j) 
k )i 


dJp \- k \ tjj) 


dn 


(j) 


dn 


(j) 


(P.„op 




<^Pr 


[fc+nd-1] / (j) 


V 

d-i Pu 


(pS: 


cF P\ 


[fc+nd-1] / (j) 


chn 


(?) 


(P u 


, 0 ), 


O') 

G M rfe ' ndXI/ i 


(V, + 00 Mfc+nd] / (j)\ 

0[fc,nd]Ji —I ~ (i) [Pi,l),- 

v an ' 




<9n 


jn)) e M l[fe+nd]|x P , 


in terms of which we write the homogenous sample matrices 


( S nd) ( i) , (Ef ) U) ) G 


md\ 0) 


nd\ O') A 
Pj 

(S[Md])^ » ( S m)^) e M I[M1|X ^ 

which allow us to define the multivariate Wronski type sample matrices 

F'f : = ((E£ d ) (0) ,,..., (E” d )( d-1 )) g R rfc ’" dXrfc ’" d , 

£[fc,nd] :=((S M ) (0) , • • •, (S M ) {d - 1} ) G 

Definition 4.2. We say that Af ki nd is a poised set if the sample matrix is non singular 

det F f ^ 0. 

Theorem 4.1. For a poised set of nodes Af k , n d C Z(Q) in the algebraic hypersurface of the prime 
polynomial 1Z the transformed orthogonal polynomials can be expressed in terms of the original ones 
as according to the quasi-determinantal expression 


p j „0) 


Pj ( j) 

UK 



( 

P \k] (*) 

\ 

TD )[fc],[fc+nd] ^ 

TP ' k ' [x) = K(xY °- 

Yj n d 

P[k+nd— 1] (®) 



V ^[k,nd\ 

-P[fc+nd] («®) 

/ 
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Proof. Proposition 4.2 gives 

^[k\,[k+nd\ ^[k,nd\ 4“ j ' ' ' j ^[k],[k-\-nd— 1]) 0 

so that 

(^[fc],[fe]) • • • j ^[k\,[k+nd— 1 ]) = — (77(A) ) [k]\k+nd}*-‘\ k ’ nd \ ) ’ 

and 7 Z(x) d TP(x) = c oP(x) gives the result. 

To discuss the existence of poised sets we allow the nodes to be complex. 
Definition 4.3. We introduce the partial derived Vandermonde matrices 


□ 


TT (p M) VlPl (p u) ,\ eC 


fc+mf— 1 ^ 


0) 


’ <9n- 

for j G {0,..., d — 1} and i G {1,..., pj }, the derived Vandermonde matrix is 

(Vf) w : = ((Vp)»,.... (Vp)«) £ C"*+" 

and t/ie multivariant confluent Vandermonde matrix 

Vf d := ((Vfc d ) (0) , (V fc nd ) (1) ..., (Vfc d ) (d_1) ) G C* fc +" 

As in the previous analysis we have E™ = S™V™, where S™ is given in (3.1), and Ker Vf d C 
KerE£ d . For Mk.nd to be poised we must requ est to V £ d to be a full colum n rank matrix; i.e., 
dim Ini Vf d = rp.., u i- Remarkably, Theorems 3.1 and 3.2 and Propositions 3.6 and 3.7 hold true for 


Vnd-lXYflA j) 


\-nd— 1 X-Vk nd 


our polynomial Q = lZ d and the corresponding confluent Vandermonde matrix Vjf, m = nd. 

Theorem 4.2. The node set J\fk,nd C C D C Z(Q) is poised if it does not exist a polynomial V lZ d , 

’ (i) / Qiy \ 

deg V < k + nd — 1, such that A ffl C Z (- r— j, for j G {0, ..., d — 1} and i G {1, ..., Pj}. 

^ dnf ' 

Proof. A vector v = (uj) i= fc 1 + " ,d_1 G (ImV^) 1 if for the corresponding polynomial V = vix ai the 


d j V o d j (x a 1Z d ) 

polynomials- 77 - cancel at Ay . Remarkably, - -r. — (p) = 0, 


1=1 


dnP 1 dnP 

and p G Z(JZ). Hence, we conclude that 

(R d )k+n d -i C (ImVf 


= 0, j = 1, ■.. ,d-l for i G {1,... ,pj} 


and, as clim(77 d )fc +ri rf_i = Nf.-i, we have full column rankness of the confluent Vandermonde matrix 
if ImV^ d = {'R' d )k+nd-V □ 

Corollary 4.1. For a poised set 

• we can not take the vectors np G M^ll such that for a given p G {1,..., d} the polynomial 
d j (K^ 


d n 


(j) 


cancels at Z(fJZ). 


we can not pick up the nodes from an algebraic hypersurface of degree less than or equal to 

LvrJ + "fl 

the following upper bounds must hold 


( 0 ) ^ 

^0 — ^k-\-n(d— l),n? 


(j) 


where dp := deg ■ 


dnP 


2 Here [yj is the floor function and gives the greatest integer less than or equal to x. 
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dHlZ p ) 

Proof. When -.cancels at Z( 1Z) then dim(Im V™) L > W+(d- P )n and the set is not poised. Given 

dnf’ 

a polynomial W, deg W < |_^Aj + n, of the described type we see that V = W d , deg V < k — 1 + m, 
is a polynomial such that cancels at Z(W) and again full column rankness is not achievable. 

All the columns in the Vandermonde block (V£ d ) ( ' 0 ' ) (we have remove the subindex because for j — 0 
there is only one and no need to distiguish among several of them) imply no directional partial 
derivatives, so that (JZ)f +nd -i D Im(Vj[ n ^)^ and the maximum achievable rank for this block is 
N k+ nd-i ~ N k+n ( d _ 1)_1 = r k+n (d-i),n- For j = 1,..., cl - 1 the columns used in the construction 

D 


of the block (V™ i )p' ) imply directional partial derivatives an d consequently 


dri- k-\-nd —1 


Im (Orf) j; hence, the maximum rank is N k+ndi - 1 — ^ k+nd _ 1 _d (S> ' 


□ 


4.2. The general case. We now consider the general situation of a polynomial in several variables, 
i.e., Q = lZ dl ■ ■ ■ 1Z d jf' where 7 Zi, deglZi = mi, are different prime polynomials; we have for the degree 
of the polynomial deg Q = m = rqdi + • ■ • + n/v(£v- As with the study of the product of N different 


prime polynomials developed in §. 


N 


we have Z(Q) — (J Z(7Zi). 


i= 1 


Definition 4.4. Consider the splitting r k , nidl+ ... +nNdN = N k+nldl+ ... +nNdN _ 1 - N k _i = \[k]\ 4-b 

N d a — lPa,j , .. 

| [k + dnidi + • • • + n]yd]y — 1] | = q > an d f or each j G {0, 1,..., d a — 1} the following 

a=l j =0 2—1 

set of different nodes 


:={ P <y>}£'' c 




d a - 1 Pi 

V := U UW’ 

3 =0 *=1 


where for a fixed a G {!,..., N} we allow for non empty intersections among sets with different values 


N 


of j; denote its union by M k ,n 1 d 1 +-+n N d N = U -^4 and pick a set of linearly independent vectors 


a= 1 


C M^'ll = (R D ) ej , pj < \[j]\, j G {0,... , d a — 1}. The associated i-th homogeneous blocks 
of the sample matrices are 


( 


(E 


nidiH- hn N d N \(a,j) 


(PP- 


dn ) 


- [k] / (a,j)\ 

AyKr) 


P[ k +mdi -1-hnjvrfjv—1] / (aj)x 

\Pi, 1 t 


(£[*,■ 


nidiH-hnjvdjv] )i 


(®j) ._ 


V 

cPP- 


dn] 


(“j) 


[fc+rtldiH- \-n N d N ] / (a,j) 


An 


(<*j) 


(Pirr), 


A J / 3 - 


M_ / (a,j) X 


An, 


\ 


(aj) 


- \-n N d N -l] , ( a,j) \ 

din {a ’ j) [P iA a,j) J 


/ 


d j Pi 


[fc+nidiH- \-n N d N ] , ( a,j ) 


An 


(®.j) 




With (Y 1 ^+-+ n ^)\ a d) G M r fc ,„ ldl+ ... +rlJv ^x^“’ j) an rf(S [fcinidl+ ... +njvdiv] )5 aJ) G Kl[fc+nidi+-+n^^]|x^ a ’ 3) j 
the homogenous sample matrices are 


-vnidiH- ( a ij) . / /yinpdiH-( a ’*?) 


- 


\ \ a ij) / 

Jl 

IV 1 \( a o) . _//vi \( a J) (yi \( a J)\ _ ]n)|[fc+nidiH-hnjvrfjvll x V 

[fc,nidiH-hnjvrfjvjJ — V V ^ rfi-f- \-n N d N ])± > ■ ■ ■ i f z - , [fe,nidi4- \-n N d N ]) pa ) t JK. 


nidiH-h^jv^ivA (°' 0 ) 


^ - fnjvdjv xE i= N 


p a,j (a,j 
i=l u i 
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and the partial multivariate Wrohski type sample matrices are 

/yinidiH- \-nffdff\ . _//ytnidi-t-hnjydjv\(a,0) /v-inidiS-hriNdjv\(a,d a —l)) 

K^k > a — > !!■■■) K^k > 

(S[fc,nidiH-fnjv<ijv])“ - = ((^[fe,nidi-+- \-n N d N ])^ ’ \ • ■ ■ j (^[fc,riidi4 - fnjvdjv])^ ’ ^)j 

v \r^d a — 1 ^~\Pa,j (ci,j) ir* , j . . 1 — 1 ^-\Pa,j ( a ij) 

that are matrices inM. rk ’ nidl+ ' +nNdN ^i=° ^ i=1 Vi and in IP +ni lH hnjv ,iVjl r= o l^i=i v i Respec¬ 
tively; finally, consider the complete sample matrices collecting all nodes for different a € {1,, N} 


-\7i\d\-\ - \-ri]\[d]\[ , 


:=((S 


nidi-I- \-n N d N 


)i, 


nidH-hnjvdjv 


)n), 


- l-n N d N ] ■ ((^-'[fc,mdi4-fnjvdjv])lj • • • j (^-‘[fc,nidi4-hnjv^jv])Jv) • 

We use the word partial in the sense that they are linked to one of the involved prime polynomials. 
We now proceed as we have done in previous situations just changing nodes and sample matrices as 
we have indicated. Then, 

Definition 4.5. We say that Afk,n 1 d 1 +-+n N d N is a poised set if the sample matrix is non singular 

7^0. 


det —tn/v<tv 


N 


Theorem 4.3. For a poised set of nodes fifk,n 1 d 1 +--+n N d N in the algebraic hypersurface |J Z(Q a ) the 

a= 1 

transformed orthogonal polynomials can be expressed in terms of the original ones as the following 
last quasi-determinantal expression 


rpr> ( ria-l( , ^' a (-^-)) ) [fc],[fc+mdi4- \-n N d N ] n 

TP[ ‘ i(!e) = e * 

( 

P [k}(x) ^ 

Y^n,idi~\ - \-n]\fdN 


-P[fc+nicM-hnjvdjv-1] (• K ) 


V -bnjvdjv] 

-P[fc+nidiH-bnjvdjv] I 3 ') / 


Proof. Proposition |4.2| gives 

^[k],[k+mdi-\ -hnjvdjv]^[fc,nidiH fnjvdjv] d~ > ^[fc],[fc-|-nidi4 bnjvdjv —1]) 0, 

so that 

(^[A;],[fc]) • • • ; ^[k},\k+n\d\ + —hnjyrfjv —1]) 

N 

= - (n(^(A))“-)) |t] , |l+ni , 1+ ... +w ,„ | SM,* + - 1 


a =1 


and lZi(x) dl ■ ■ ■ lZ N (x) dN TP(x ) = uP(x) gives the result. 


□ 


As in previous discussion we shift from the field M to its algebraic closure C, laying the algebraic 
hypersurfaces and nodes in the D-dimensional complex affine space. 


Definition 4.6. We introduce the partial derived Vandermonde matrices 

-( 


, ■ \ / Ai-yA+’U'AH-1 -n N d N ) , f)j v [fc+mdi4 - \-n N d N ] 

f\?nidi-|- \-n N d N \(a,j) ,_( L V u A (^( a O) W 

)i — \ „ ( a : j ) yPi, 1 „ (a.i) ^ ) ] , 


dn 


dn 


(®d) 


.(o.i) 


that belong to C Nk+n ^ + "' +n N d N - lX V '' . For j e {0,..., d a — 1} and i G {1,..., p a ,j}, the derive 0 d 
Vandermonde matrix is 


(ynidH-b n N d N ^(a,j) ._^ynidi-+ b n N d N ^(ji,j) -tnjvdjv^aj)) ^ ^A fc+ „ 1( j 1+ ... + „ ivrdjv ._ix2i=i iA 

AT \^ da , — 1 ( 

and the partial multivariant confluent Vandermonde matrix in C fc+r *i d i+ ■ +n Jv d iv- 1 Ai=o ^»=i ^ 

^ynidi-f-bnjvd/v) ._ A‘y n i'tiA-bnjvdjv^(a,0) ^ynidiH-bnjvd/v^(a,l) ^ynidi-f-bnjvdjv)(a,d a —1)\ 
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Finally, the complete confluent Vandermonde matrix in G (C Nk + n i d i+-+^N d N - lXrk ^i d i+-+n N d N ^ 

._ -friivdjv^ (ymdi4-h n N d N \j 


-p«ldl4- \-n N d N 


Again we find the factorization x£i*+-+"^ = s nid 1 +-+n N d NV n 1 d 1 +-+n N d N ^ with gm as in 
so that KerV” ldl+ - +njvdjv C Ker ^fl dl+ - +nNdN . For Af^+-+^d N to be poised we must request 
to Yfl dl + ' +n N d N |- 0 j iaye f u jj co i umn ran k matrix: dimIm V k ldl+ '-+n N d N _ r k , nidl+ ... +nNdN . Again, 
Theorems 3T and A2 and Propositions AG and A7 hold true for our polynomial Q = lZ d and the 
corresponding confluent Vandermonde matrix m = nidi + • • • + n^dAr. 

Theorem 4.4. The node set Mk,n 1 d 1 +-+n N d N C C D is poised if it does not exist a polynomial V 


ria=i deg V < k + nidi + • • • + njydN — 1, such that J\f- 


Oh) 


C 


z( 


HV 


' dn. 


Oh) 


; for a = 1,..., N, 


j e {1,... ,d a - 1} and i G {1,... ,p a j}- 


N, 


Proof. A vector v = {vf) i= \ 


k+nidi~\ - \-n N d N — 1 


6 (Im V t 
d’v 


nidi 4-hnjvdjvlT 


if for the corresponding polynomial 


V = J2i^i’ l i dl + '" + " Nd N 1 viX a ‘ the polynomials ° cancel at J\fj a " j> . Notice that the polynomial 


V = x a nf-i T^h b i s such . ., 

6-1 6 dnf- j) 

{1,..., d a — 1} and i G {1,..., p a ,j}- Hence, 

N 

U n 


d j V 


dn 

do cancel at \jfl =l Z{JZb) for j — 0 and also at Z(lZ a ) for j G 


a= 1 


k+n\di~\ - \-nppdN — 1 


C (Im Vi 


nidiH- \-rt]yd N \± 


)\ 


which considered at the light of the condition dim(J|^ 1 lZ da )k+ ni d 1 +-+n N d N - 1 = N k -i implies that 
no further constraint can be allowed or the full column rank would not be achievable. Thus, the set 


is poised if Im V k 


nidi 4-(-njvdjv 


= {iz a 


'k+nid\-\ - \-npfdjsf— 1' 


Proposition 4.3. In order to have a poised set the polynomial 


dflTZl 


pi 


dn 


(®h) 


□ 


KW) 

- can not cancel at 


Z{1Z\ ■ ■ ■ TZn) for 0 < pi < di, - - - , 0 < pn < dx. Moreover, when the set of nodes is poised we can 
ensure the following bounds for the node subset cardinals 

k 


(4.1) 

(4.2) 

(4.3) 
Here 


\K\ > 


d n 


+ n a , 


\N {a ’° A ^fc4-nidi4-hnjvdjv— n a ,n a j 

I A/'(“h)| < r . 

I ® I — fc4-nidi4- \-nffdff—dp’ J ,n 


d(°h) ;= deg ^(^a Q ) 


dn) 


(®h) 


and the function [x] gives the smallest integer > x. 


d'HV'F ■ ■ ■ V PN ) 

Proof. If-- —ttt— cancels at Z(TV) then 

dn [ p 

dim ((ImK 1 *^ 4 "'") 1 ) > Y + „ 1 * + ... + „rf„- (n „ 1+ ... +w *) 
and the set is not poised. 
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To prove (4.1) let us consider first order polynomials 7 t 4 -“’^(z) which cancels at p\i, To (Pu^) = 
0, and construct the polynomial hl a = ( j=\ ,. tda 7 Vi), degll a = |A/" a |. Then, the polynomial 


V = IZf 1 ■ ■ ■ ^a-i'^'a+i ''' deg V = nidi H-f ryyh/v + (|-A/" a | — n a )d a , has the nodes among 

r\j T T 

its zeroes J\f™' d '+'" +nNdN Z(V) and ^ do cancel at for all b G {1,..., TV}. Thus, we 

should request nidi + • • • + uncIn + (|A/" a | — n a )d a > k + nidi + • • • + n^d^ ; i.e., (|A/" a | — n a )d a > k 
and the result follows. 

For (4.2) observe that all the columns in the Vandermonde block (v” ldl+ ' +,lJvdjv )( a .o) imply no 
directional partial derivatives and are evaluated at nodes which belong to Z(lZ a ). Therefore, 


■» d a — 


1 =1 .^ 


(w.) 


k-\-n\d\ H-h^-AT^iV 


_! D Im(Vr dl+ '" +nivdiV )i 0) 


and the maximum rank achievable for this block is N k+nidl+ ... +riNdN _ x - N k+nidl+ ... +nNdN _ na _i = 
r k+nidl+ - + n NdN -na,n a - Tlie columns of the block (y£ ldl +"'+ n iv-djv)(“j) are vinculated to 

and consequently (^J 1 )i +fe+ „ ldl+ ... +njvdjv _ 1 2 (l m (V^ +nidl+ '" +njvdjv )? ) ); hence, the maximum 
possible rank for this block is N k+k+nidl+ ... +nNdN _i - . □ 
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